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ABSTRACT 

We analyze a (120 Mpc) 3 adaptive mesh refinement hydrodynamic simulation 
that contains a higher-resolution 31 x 31 x 35/i -3 Mpc subvolume centered on a ~30 
Mpc diameter void. Our detailed ~1 kpc resolution allows us to identify 1300 galaxies 
within this void to a limiting halo mass of ~ 10 10 A/q. Nearly 1000 galaxies are found to 
be in underdense regions, with 300 galaxies residing in regions less than half the mean 
density of the simulation volume. We construct mock observations of the stellar and gas 
properties of these systems, and reproduce the range of colors and luminosities observed 
in the SDSS for nearby {z < 0.03) galaxies. We find no trends with density for the most 
luminous (M r < —18) galaxies, however our dwarf void galaxies (M r > —16), though 
they are less reliably resolved, typically appear bluer, with higher rates of star formation 
and specific star formation and lower mean stellar ages than galaxies in average density 
environments. We find a significant population of low luminosity (M r ~ —14) dwarf 
galaxies that is preferentially located in low density regions and specifically in the void 
center. This population may help to reduce, but not remove, the discrepancy between 
the predicted and observed number of void galaxies. 



Subject headings: galaxies: evolution — cosmology: theory 
scale structure of universe — methods: numerical 



hydrodynamics — large- 



Introduction 



ACDM has proven fairly robus t when compared with a wide range of existing observation s 
(|Lacev fc Colelll993l : ICen et al.l[l~99l : Izhang et al.1ll995l : ISpringel et aUbood : komatsu etldlboilh . 
Simulations can reproduce the galaxy clustering properties we observe in large redshift surveys, as 
well as the size of galaxy cluster s. In addition, ACMD predicts hiera rchical growth as the dominant 



mechanism in galaxy evolution ([Peebles 



1969 



White k Frenklll99ll 



However, discrepancies with observation remain, particularly on individual galaxy scales. The- 
ory predicts a cuspy core dominate the shape of the dark matter density profile, while observations 
show that a wide range of halo sizes, from galaxy clusters to dwarf spheroidals, instead have cored 
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profiles (jTvson et al.lll998l ; iKlevna et al.ll2003l ). Simulations have difficulty reproducing exponen- 
tial, bulgeless disks, and observations w ithin the Local Volume disagree with predictio ns of where 
we should find the largest disk galaxies ([Peebles &; Nusserll2010l ; iKormendv et al.ll2010l ). The miss- 
ing satellite problem de scribes from the overprediction of sma l l dark matter halos c lustering around 
the Milky Way Galaxy (jKlypin et al.lll999l : iMoore et al.lfl999l : Isimon k Gehall2007n . The void phe- 
nomenon refers to the ov erprediction o f the number of low-mass halos existing within voids in 



cosmological simulations ([Peebles 



2001 



One key to solving these problems may come from the 
inclusion of sufficient physics to properly simulate the baryon content of galaxies, which does not 
dominate their mass but contributes significantly to small scale dynamics and observables, such as 
stars or gas. 

Early cosmological simulations of voids were based on dark matter only, N-body simula- 
tions (e.gjRvden fc Turner! 1 19841: Iwhite et all 1 19871 : Little fc Weinberg! 1 1994 : IVogelev et all 1 1994* : 
Mathis fc Whitdbooj ). These authors compared the dark matter distribution with the void statis- 



tics in large surveys such as the CfA redshift catalogs and the Las Campanas redshift survey, in 
terms of their sizes and abundan ces. Using semi-a nalytic prescriptions for galaxy formation on 
outputs from N-body simulations. iPatiri et al.l (|2006l ) studied the colors and specific star formation 
rates of void galaxies, finding that there are more blue galaxi es in void but no system atic differences 
between void galaxies and the general galaxy population. ICeccarelli et al.l (|2006l ) used a similar 
method to study void dynamics and the effect of redshift distortions in void i dentification. Th e 
only cosmological hydrodynamic simulations ru n to date focusing on voids are IViel et al.l (|2008l ). 
who analyzed the void statistics at z ~ 2, and iHoeft et al.1 (|2006l ) . who investigated the effect of 
the cosmological UV background on the formation of dwarf galaxies in voids. 

We have selected a void region from within a full (120 h~ l Mpc) 3 cosmological simulation, and 
we examine with moderate resolution the dark matter, gas and stars of a large sample of galaxies 
located within and around this large void. We employ an adaptive mesh refinement (AMR) code 
to reproduce the gas physics on a scale of 1 kpc. Observationally, void galaxies are distinguishable 
from galaxies in average or overdense environments as they are typically less luminous, and at fixed 



luminosity are bluer, with higher rates of star formation and specific star formation (IRojas et al 



2004 . 120051 ). This is not strongly reflected in their gas content, though they are gen e rally gas-rich. 



as th eir total mass in H I is fairly typical for their luminosities (jSzomoru et al.lll996l : iKreckel et al 

bonh . 

We compare the integrated properties of void galaxies from our simulation, described in Section 
[2j both with observed void galaxies and galaxies in higher density regions within this simulation. 
We present our results in Section |4] and our conclusions are summarized in Section [5j 
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2. Simulation Initial Conditions and Physical Processes 



We perform co s mological simulations using the adaptive me sh refinement Eularian hydro code 



Enzo (jBrvanl 1 1993 ; iNorman &; Bryan 



1999 



O'Shea et al.l 12004 ). with a periodic box of size 120 



h 1 Mpc comoving on a side and cosmological parameters taken from the WMAP5 ACDM results 
combined wit h measurements of T ype la supernovae and Baryon Acoustic Oscillations in the galaxy 
distribution (|ffinshaw et al.1 l2009h : A , ^6, h, a 8 , n s ) = (0.279, 0.721, 0.0462, 0.701, 0.817, 

0.960). We first run a low-resolution simulation with a uniform dark matter particle mass of 7.6 
xlO 10 Mq, 128 3 root grid cells, and only 4 levels of refinement from an initial redshift of z = 99 to 
0. Based on the z = output of th i s simu lation, a large void region having a diameter of ~30 h^ 1 
Mpc was identified, (see also ICenl (|2010l )). Tracing the dark matter particles back to the initial 
redshift of z = 99 showed that the void region expands with time in terms of comoving volume, 
as expected. To achieve high mass and spatial resolution, we then use the multimass initialization 
technique and employ three nested volumes with successive particle masses decreased by a factor 
of 8. Hence, the innermost 31 x 31 x 35/i~ 3 Mpc comoving volume, sufficiently large to contain 
the identified void volume at z = 0, has a dark matter particle mass of 1.5 xlO 8 M Q . Within 
this innermost nested volume, hydrodynamic refinements were allowed beginning with a root grid 
cell size of 937 h~ l kpc and a maximum refinement level of / max = 10, resulting in a maximum 
resolution of 0.916 h^ 1 kpc at z = 0. 



The simulation includes a metagalactic U V background (IHaardt fc Madau 



of photoelectric and photoionization h eating (I Abbott 



hydrogen (jAbe 



et al 



of UV rad iation by neutral hydrogen (|Cen et al.ll2005l ). They also includ e cooling due to molecular 



19821 ; IJoung Mac Low 



1996), a diffuse form 



20061 ). and shielding 



19971") and metallicit y-dependent radiative cooling (jCen et al.lll995l ) extended 



down to 10 K (IDalgarno fc McCravlll972l ). 



Inside the innermost nested volume where AMR is allowed, additional physical processes of 
star formation and feedback are imple mented. Star part i cles a re created in cells that satisfy a set of 
criteria for star formation proposed bv lCen <fe Ostrikerl (|1992l ). which requires that the gas within 
that cell be contracting, cooling rapidly and gravitationally unstable. Under these conditions, a 
stellar particle of mass m* = c*m gas Ai/i* is created to replace gas from that cell, and is tagged with 
its initial mass, creation time, and metallicity. Here, At is the time step, t* = max(tdyn> 3 x 10 6 yr), 
t dyn = y/3~n '/ (32Gptot) is the dynamical time of the cell, m gas is the baryonic gas mass in the 
cell, and c* = 0.03 is the star formation efficiency. Star particles typically have an ini t ial m ass of 
~ 5 x 1O 6 M0. Star formation and supernovae feedback are modeled following Cen et al. (j2005l ) with 
a supernovae efficiency of esN = 10~ 5 . Star particle masses decay slightly as feedback energy and 
ejected metals are distributed into the 27 local gas cells, weighted by the specific volume of each cell, 
and centered at the star particle in question. The temporal release of metal enriched gas and thermal 
energy at time t has the following form: f(t,ti,t dyn ) = (l/t dyn )[(t - U)/t dyn ] exp[-(i - ti)/t dyn ], 
where tj is the formation time of a given star particle. The metal enrichment inside gala xies and in 



the in tergalatic medium (IG M) is followed self-consistently in a spatially resolved fashion (|Cen et al 



20051 ) . Recently, ICenl (|2010l ) used a simulation of the same system with a higher spatial resolution 
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(by a factor of 2) to study the nature of damped Ly-a systems. 



3. Mock Observations within the Simulation Volume 



We identify yirializ ed objects in our high-resolution simulations using the HOP algorithm 
(jEisenstein &: Hutl Il998l ) with the threshold parameter (5 uter) of 125. For identifying galaxies, 
we used higher <5 ou t er values of 10 3 and 10 4 to additionally find subhalos located within virialized 
objects. Galaxies are excluded that are outside of the central, high-resolution ~ (30 /i" 1 Mpc) 3 
region, or that contain any coarse dark matter particles more massive than the highest particle 
resolution achieved. We also exclude any dark matter halos with fewer than 100 particles within 
the virialized region to ensure sufficient resolution of the physical balance between pressure and 
self-gravity within the simulated galaxies. For this work, we have extracted the hydrodynamic data 
around every identified galaxy at a uniform resolution corresponding to level 8, or a physical scale 
of 3.66 hr x kpc at z = 0, as it is sufficient for our examination of the integrated properties in these 
systems and significantly speeds our analysis. 

The light distribution is computed from the star particles usi ng the Galaxy Isochrone Syn- 
thesis Spectral Evolution Library (GISSEL) stellar synthesis code (jBruzual &: Charlotll2003l ). We 
calculate the luminos ities of the simulated galaxies in the five Sloan Digital Sky Survey (SDSS, 
Abazajian et al.ll2009l ) bands (ugriz). Internal extinction corrections were applied by considering 
the mass column density of metals, T,z, along the line of sight within the virial radius for a ran- 
domly chosen viewing direction along the simulation volu me axes. Starting with the observational 
relation determined in the Milky Way at solar metallicity (IBinney &: Merrifieldlll998l ). 



Ax 



N(H t , 



1.9 x 10 21 cm" 1 



mag, 



(1) 



we allow for s ome dependence on metallicity, and scale for the fraction of refractory elements, fp e 
(jviadilolbooi ). Thus we find that 

/Fe 



F m p 4 x 10 19 cm" 2 



mag. 



(2) 



We choose a scaling factor, F = 1.5, to match our simulated extinction with the observational 
relation (Equation [T]) for those simulated galaxies with solar metallicity. We also considered the 
observational dust to gas relation determined in the Large Magellanic Cloud, which has a muc h 
larger uncertainty, as compared to 1/3-1/2 solar metallicity simulated galaxies (iKoornneefl Il982l ) . 
We find a best fit with a factor of F ~ 3, however a factor of F = 1.5 is still in agreement to 
within the errors, which are quite large, and we adopt this value for the remainder of this paper. 
Extinctions for each star particle w e re th en scaled to the central wavelength for each of the five 
SDSS bands following ICalzetti et al.l (|2000l ). We define the integrated stellar parameters (i.e. mass, 
luminosity, star formation rate) of each simulated galaxy to be the sum of the properties of the 
star particles located within 15% of the virial radius of the galaxy at a given redshift. 
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The dark matter density field was determined by calculating the mean dark matter density 
for the entire simulation volume, gridding all dark matter particles to a 1 h~ l Mpc grid, then 
applying a 5 h~ l Mpc three-dimensional Gaussian smoothing filter. The smoothing length was 



chosen to agree with t he galaxy correlation length observed in large redshift surveys (|Jing et al 



19981 ; iPark et al.l 120071 ) . We also note that our choice of smoothing length assigns a reasonable 
dark matter density contrast of 5p/p < —0.5 to all galaxies where the distance to the third nearest 
neighbo r is greater than 7 h~ x Mpc, which was the more restrictive void galaxy selection technique 
used by iRoias et al.l ((200J) . 



4. Results 

Of a total of 1268 galaxies with high resolution hydrodynamic simulation, none are found at 
positions with a dark matter density contrast (see Section [3]) of more than 1 as we include only 
filaments and walls bounding the void but no massive clusters. 950 galaxies are found at locations 
with Sp/p < 0, from which we form a void sample (VS) of 648 galaxies with —0.5 < Sp/p < 0, and 
302 that are in regions with 5p/p < —0.5 form a low-density void sample (LVS). There are also 
318 galaxies that make up a non-void sample (NV) at roughly average densities (0 < Sp/p < 1) 
on the void boundaries (Table HJ Figur e CD) . We make these di visions in part to facilitate direct 



comparison with observed void samples (jRojas et al.ll2004l . 120051 ). We can also define the distance 



R, from the void center at [64, 71, 45] Mpc, as identified within the full 120 h" 1 Mpc box, and 
find that only two galaxies resolved by the simulation are within 5.5 h^ 1 Mpc and no NV galaxies 
are within 18 h~ l Mpc. We use this to additionally define a void center sample (VC) of 176 galaxies 
within 18 h~ l Mpc. Simulated observations were made with a typical H I column density detection 
limit of 1 x 10 19 cm -2 , and resulting integrated H I column density contours overlaid on simulated 
observations of the stellar luminosity (Figure [2]) . 

We apply a Kolmogorov-Smirnov test to determine if the luminosity distribution of these 
samples could be drawn from the same parent population. The probability that the LVS or VC 
samples are drawn from the same population as the NV sample is fairly low (P < 0.03), however 
there is a very strong probability (P = 0.847791) that the VS sample is drawn from the same 
population as the NV sample. This suggests that there is a sharp distinction between the galaxy 
populations in the deepest underdensities and those in more moderately underdense regions. This 
is most apparent in the large population of low luminosity galaxies found preferentially in the void 
center (see Section l4~5j) . 



4.1. Luminosity Function 



As the most luminous galaxies have been more robustly simulated throughout their merger 
history, we expect that this population will most closely match observations. Indeed, the observed 
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void galaxy luminosity function is very well reproduce d (Figure [3|). Here we have divided our sample 
into two at Sp/p < —0.4, following iHovle et al,l (|2005l ). and contrast the lower density void galaxies 
with the higher density 'wall' galaxies. Our simulated void sample very closely matches their 
observationally determined Schechter function fit. Our 'wall' sample contains too few galaxies, 
however our volume excludes any significant high density regions and a significant number of 
galaxies contaminated by coarse particles. The turnover in both samples at M r = —18.5 indicates 
the completeness limits in our simulated galaxy samples, and we expect that increased resolution 
would not significantly change our results for galaxies above this luminosity threshold. 



4.2. Colors 



Compared to observations, our simulated void galaxy colors appear consistent, however they do 
not differ greatly from the non-void gal axies. Table [ 2] com pares the mean properties of our galaxy 
samples with the photometric study by iRoias et al.l (|2004h , who construct a void galaxy sample of 
those galaxies with a third nearest neighbor distance greater than 7 hr 1 Mpc and a wall sample of 
the remaining galaxies. Duplicating the magnitude limit (—17.77 > M r > —22.5) for their distant 
galaxy sample, and subsequently dividing it into a bright (M r < —19.5) and faint (M r > —19.5) 
subsamples, we see no consistent trend with density in our galaxy colors, but a strong dependence 
on magnitude. I n all magnitude limited s amples, we note that the NV galaxies are too blue when 
compared to the IRoias et al.l (|2004l . 120051 ) wall sample. 



In addition, none of the samples exhibit bimodality of galaxy co lors with a red sequence and 
blue cloud (Figured]) such as is observed in both void and field galaxies (Ivon Benda-Beckmann fc Miiller 
20081 ). We compare the LVS with the geometrically sele cted Void Gal a xy Su rvey (VGS, Kreckel et 
al. 2011, in preparation; selection criteria described in iKreckel et al.l (|201ll )). and the NV with a 
magnitude- limited sample of SDSS galaxies selected at redshifts 0.1 < z < 0.3. Both comparisons 
show general agreement between simulation colors and observations down to the observational limit 
M r ~ —16. The choice of a lower scaling factor, F, in Equation [2] corresponds to more extinction 
and redder galaxies, however it affects void and non-void colors similarly, does not affect the bi- 
modality, and produces an excess of low-luminosity red galaxies compared to the VGS. The lack 
of bimodality in color in the NV sample (Figure [5]) suggests that the homogeneous distribution of 
galaxy colors is independent of density. 

While the absolute value calculated for our galaxy colors may not be entirely reliable, we may 
still co nsider relative differences between the void and non-void samples. Ivon Benda-Beckmann &: Miiller 
(|2008l ) detect a slight blueward shift of the blue cloud galaxies that are deeper inside the void at 
fixed luminosity. Considering the mean color at fixed luminosity (Figure [5j left) we see no such 
shift at the high-end luminosities they were considering (—19.2 > Mb > —20). They have fit their 
color distribution with two gaussia ns, which we cannot replicate as we see no bimodality, however 

we still find no trend with density (Figure 



2.2 rtStrateva et al. 



2001 



dividing the sample at u — \ 
El right). We do see indications of a blue shift at fainter magnitudes, M r > —18. This is presumably 
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due to the higher rates of star formation in galaxies at these densities, described below. 



4.3. H I Properties 



The H I distribution typically forms a disk, however the ~5 kpc resolution used for the extrac- 
tion of the hydro variables is insufficient to attempt any analysis of the internal gas disk kinematics. 
The H I surface density profiles we find rise too steeply at the outskirts of the disks, often achiev- 
ing fairly high colu mn densities of ~ 10 21 cm -2 far outside the stellar disk (Figure [21 right) (cf 



Swaters et al.1 12002 ). This inaccurate density profile persists even if we are careful to include in 



the integrated H I map only that H I that is above a detectable level (1 x 10 cm ) for modern 
instruments. We believe our measured H I masses to be an upper limit based on two limitations 
with this simulation. First, the conversion of H I to H2 through formation of H2 on dust grains was 
not properly modeled even at this resolution. This should result in an overestimate of the H 1 mass 
in the highest density, central regions of the gas disk. Second, interstellar UV ionizing radiation 
by individual galaxies was neglected, which would have reduced H 1 in the outskirts of the disk, 
because it is likely that local ionizing radiation is larger than the meta-galactic background. 

The resulting overestimate of the H 1 mass is apparent in comparing the H 1 mass to light 
ratio to observations (Figure [6l left). It shows the same decreasing trend with increasing luminos- 
i ty, however the overall value is about a facto r of three too high c ompared with typic al galaxies 



(jVerheiien Sancisi 



2001 



Swaters et al.ll2002l ) or with void galaxies (jKreckel et al.ll201ll ) . We note 



that there is a slight trend with density among our three density-selected samples (Figure [61 right), 
however when we examine the specifics of the trend, particularly for low luminosity and dwarf 
galaxies (Figure [7}, i t is not very pronounced and certainly much weaker than the trend observed 
for dwarf galaxies by iHuchtmeier et al.l (j 19971 ) . We see no trend as a function of distance from the 
void center for our VC galaxies. 



4.4. Star Formation and Stellar Properties 

Although there is a no strong change in the star formation rate (SFR) at lower densities (Table 
[2]), it is clear from Figure [8] that there is a strong dependence of SFR on luminosity. Thus we also 
consider the SFR scaled to the stellar mass (specific SFR, S-SFR), which is significantly increased 
for observed void galaxies, but unchanged within the error in the mean for our simulation (Table 
[2]). We do, however, see significantly higher SFR and S-SFR for those lowest luminosity dwarf 
galaxies with M r > —16 in underdense regions when compared to our NV galaxies (Figure [8] and 
Figure left). The si mulation does not re produce the observed trend for increased SFR per H I 



mass at lower densities (jKreckel et al.ll201ll ). and there is good agreement between densities for all 
but the faintest dwarf galaxies (Figure [9|) , however our H 1 masses are only upper limits and our 
sample is not complete for low mass galaxies. 



- 8 - 



We find no trend in the mass-metallicity relation with density (Figure [T0|) , or with distance 
from the void center. This is likely because metal enrichment is mostly due to self-enrichment, 
and thus depends more on halo mass than large s cale environment. We reproduce the observed 
increase in metallic ity for more massive galaxies ([Lequeux et al,l 11979 ; iGarnett fc Shields! 119871 ; 
Zaritskv et al.l Il994l ). We also reproduce the general shape in the mass-metalli city relation as 



calcu lated using the gas-phase oxygen abundance measurement of SDSS galaxies (jTremonti et al 



2004 ). however their survey measures only the 3" of the observed galaxies, ^25% of the total light , 



in the central region where galaxy metallicity is typically higher (|Searldll97ll ; IVilchez et al.lll988l ). 
This aperture bias may result in an overestimate of the to tal metal abundance by as much as a 
factor of 2, though it should not affect the general trend (jTremonti et al.l 12004 ) . In comparing 



with t heir measured oxygen abundance, 12 + l og(Q/H), we have normali zed the iTremonti et al 



(|2004l ) relation to solar metallicity, Zq, at 8.69 (jAllende Prieto et al 



2001 



reasonably well reproduce the Milky Way, with Z at M* 
at M* = 5 x 1O 9 M . 



Our measured values 
5 x 10 10 M Q , and the LMC, with Z /2 



4.5. Excess of Low Luminosity Galaxies 



One surprising discovery is the population of low luminosity galaxies apparent below the ob- 
servational survey limit in the color-magnitude diagram (Figure 13]) • While these galaxies do inhabit 
the lowest mass dark matter halos, we believe that this excess is not a resolution effect as it persists 
into halos ~3 times more massive than we can reliably resolve and has a strong trend with density. 
Figure [TT] shows we have a slight excess of Mr ~ — 14 galaxies for our NV, a noticeable excess for 
the VS, and a substantial excess for the LVS. In the VC this excess is particularly notable. This 
excess may be due in part to the youth of galaxies in the underdense regions (Figure fl~2|) , although 
at M r = — 14 the mean ages are about the same. We note that while these galaxies are sufficiently 
resolved at z = 0, the galaxies that formed them earlier in their merger history presumably were 
not, so their integrated properties are correspondingly less reliable than in more massive galaxies. 
For example, these galaxies are generally gas rich, with masses ~ 10 9 M Q , which may be due to 
decreased gas consumption at earlier times resulting from simulation resolution limits. 

While complete volume limited samples of dwarf galaxies are not yet possible, comparison 
with observations show that our simulated dwarf galaxies are distinctly too red (Figure fl~3"j) . Here 
we compar e with dwarf galaxies within the Local Volum e that have integrated B and V band ob- 



servations ( Makarova 



1999; 



the SDSS redshift survey (jGeha et al 



Makarova et al 



2002 



200 



3T 



20091 ) . and with dwarf galaxies identified within 



We convert ugzri colors for the SDS S and simu- 



lated galaxy samples to Johnson-Cousins UBVRI colors following the prescription of I Jester et al 



(|2005l ). In Figure U] (right), when comparing our NV galaxy colors with typical SDSS galaxies, the 
simulated low luminosity galaxy population appears almost as an extension of the red sequence 
while the observed dwarf galaxies (Figure [T3]) appear more continuous with the blue cloud. While 
the observational samples are not complete in any way, there is no bias for color in their selection 
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and both samples should be sensitive to a redder dwarf galaxy population such as we simulate. 
This suggests that our simulated colors for the lowest mass galaxies, which contain very few star 
particles, are not realistic, presumably due to resolution effects earlier in their merger histories. 

Despite the discrepancy in colors, the increased number of simulated dwarf galaxies in the low- 
est density regions may present a way to reduce the discrepancy be tween the number of dark matter 
halos predicted within voids and the number of galaxies observed (|Peeblesll200ll ). iTinker fc Conrov 
(|2009l ) suggested a sufficiently tailored halo occupation distribution could solve this problem, and 
predict that the maximum halo mass (and thus luminosity) will increase as a function of distance 
from the void center. We see no such relation for VC galaxy luminosities, though we do see a slight 
trend with dark matter halo mass but not nearly as steep a cutoff (Figure [1 



Comparing the number density of void galaxies in our simulation, which likely a lower limit 
due to imperfect resolution, with the most sensitive observational surveys we still over-predict 
the number of void galaxies. If we consider observations of the Local Volum e, where the galaxy 
samp le within 8 Mpc of us is thought to be fairly complete to Mg = -12 (jTikhonov fc Klvpin 
20091 ) . and consider its overlap with the Local Void, which has a center 23 Mpc away from us and 
extends roughly 20 Mpc in radius ( Tullv et al.ll2008l ). we find one void galaxy per 394 Mpc 3 . In 
our simulation the luminosity limit is lower, but from our VC sample we expect a similar volume 
should yield about 6 galax ies. This is somewhat smaller than the factor of ten discrepancy found by 
Tikhonov &i Klypinl (|2009l ). but still quite pronounced. It is still possible that these dwarf galaxies 
in the voids are low surface brightness and thus have been largely missed by existing surveys. 
Future simulations with increased stellar resolution may be able to constrain the surface brightness 
of these systems. 

The one galaxy located in both the Local Volume and the Local Void, KK 246, has Mg = 
-13.7, a coIoilj of B — R = 1, and a significant H I mass (Kreckel et al. 2011, AJ submitted). 
Again it is too blue, but otherwise it is in good agreement with the expected properties of our low 
luminosity void galaxy population. All sky H I surveys are on the verge of reaching the sensitivity 
needed to detect these galaxies. HIPASS can detect galaxies with 1O 8 M to distances of ~10 Mpc, 
just within the edge of the Local Void. ALFALFA reported no detections within the Pisces- Perseus 
foreground void with a detection limit of 1O 8 M0 ( Saintonge et al.ll2008l ). however considering the 
discrepancy between the H I masses of our simulated galaxies compared to observations, detailed 
examination of closer voids within 20 Mpc where we are sensitive to ~ 10 7 Mq may be necessary. 



5. Conclusion 

We have run cosmological simulations using the adaptive mesh refinement Eularian hydro code, 
Enzo, with a moderate resolution refined region 31 x 31 x 35/i~ 3 Mpc centered on a large ~ 30 h 



The mean color for our simulated low luminosity galaxies with Ms > —16 is B — R ~ 1.3. 
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Mpc diameter void. We have constructed mock observations of the integrated properties (color, 
star formation rate, gas content, etc.) in a sample of roughly 900 void galaxies, 300 in especially 
low density regions and 200 located well within the void center, and 300 additional galaxies in the 
surrounding average density filaments. 

We fairly accurately reproduce the range of galaxy luminosities and colors observed, but are 
unable to reproduce the observed bimodality in galaxy colors. At low luminosities (M r > —16), 
where we more affected by resolution effect from earlier in the galaxy merger history and beyond 
the scope of most observational studies, we see that void galaxies are bluer with higher SFR and 
S-SFR, but we see no trend for the brightest galaxies (M r < —18). 

The gas content of all our simulated galaxies is anomalously high, roughly 3 times what is 
expected for their luminosity when compared with observations, and their radial extent too large, 
both of which we interpret as a result of limitations in the implementation of chemistry and radiative 
transfer in the simulations. Taking that into account, we do reproduce the trend for fainter galaxies 
to have relatively more H I for their luminosity, but find no trend with density or distance from 
the void center. We reproduce the general mass-metallicity trend observed in SDSS galaxies, and 
we see no trend with density. 

We do not see the stro ng segregation of the most luminous galaxies to the void edges predicted 



by lTinker fc Conroyl (|2009l ). however we do see that the most massive dark matter halos avoid the 
void center. We observe a significant population of low luminosity (M r ~ —14) dwarf galaxies that 
is preferentially located in low density regions and specifically in the void center. This population is 
too faint to be detected in most large scale redshift surveys, however studies of the Local Void may 
soon become sensitive to this pop ulation, where we do observe a relatively lower median luminosity 



( Nasonova fc Karachentsev 201ol ). 
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Table 1:: Categories defined for comparison between simulated galaxy populations. 



Category 


Description 


Definition 


Number of Galaxies 


NV 


non-void 


< Sp/p < 1 


318 


vs 


moderate underdensity 


-0.5 < Sp/p < 


648 


LVS 


lowest density 


5p/p < -0.5 


302 


vc 


void center 


R < 18/j- 1 Mpc 


176 



Table 2. Average color and star formation properties o f our s imul ated galaxies as com pared with 
the magnitude limited samples of iRoias et al.l (J2004J) and lRoias et al.l (|2005l ). 



Property 




LVS 


VS NV 


Rojas Void 


Rojas wall 


-17.77 >M r > -22.5 


g-r 
SFR 

S-SFR (xKT 




0.586 ± 0.0144 
0.988 ± 0.0838 
12.9 ± 1.17 


0.600 ± 0.00930 0.589 ± 0.0133 
1.80 ± 0.172 2.57 ± 0.394 
15.0 ± 1.65 14.3 ± 1.80 


0.615 ± 0.007 
0.734 ± 0.025 
3.744 ± 0.108 


0.719 ± 0.002 
0.747 ± 0.007 
2.629 ± 0.034 


-19.5 >M r > -22.5 


g -r 
SFR 

S-SFR (xlCT 


n ) 


0.593 ± 0.0234 
1.99 ± 0.318 
8.31 ± 1.02 


0.656 ± 0.0145 0.666 ± 0.0169 
4.30 ± 0.5309 6.54 ± 1.01 
8.30 ± 0.753 12.1 ± 4.09 


0.686 ± 0.009 
1.136 ± 0.063 
3.133 ± 0.169 


0.765 ± 0.002 
0.920 ± 0.016 
2.137 ± 0.004 


-17.77 >M r > -19.5 


g-r 
SFR 

S-SFR (xl(T 




0.585 ± 0.0169 
0.770 ± 0.0606 
13.9 ± 1.39 


0.579 ± 0.0113 0.550 ± 0.0171 
0.838 ± 0.0544 0.618 ± 0.0548 
17.6 ± 2.25 15.3 ± 1.79 


0.567 ± 0.009 
0.508 ± 0.0237 
4.146 ± 0.137 


0.645 ± 0.003 
0.530 ± 0.009 
3.349 ± 0.005 
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Fig. 1. — : Slice of thickness 1 h Mpc through the simulation volume smoothed dark matter 
density field. The high resolution region (white box) contains a large fraction of the large central 
void, and the positions of NV (black circles) galaxies trace the higher density regions while the LVS 
(blue circles) galaxies are located well within the void. 
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Fig. 2. — : Simulated observations of LVS galaxies, with sample observed void galaxy (bottom 
right) for reference. Greyscale indicates g-band emission, contours show H I column densities of 



cm 2 with a maximum of 1.25 x 10 21 cm 2 . The observed 



5 x 10 cm plus increments of 10 
void galaxy VGS_32 (bottom right), described in lKreckel et all (|201ll ). shows the H I beam shown 
in the lower left. The simulated and observed galaxies have roughly the same luminosity and H I 
extent, and in general the agreement is good. However, the simulated galaxies have significantly 
more low column density H I. 
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Fig. 3. — : Luminosity function of simulated void and 'wall' galaxies. Following iHoyle et alJ (|2005l ). 
we split our sample into void galaxies with Sp/p < —0.4 (triangles) and identify the remainder as 
'wall' galaxies (circles), with the respective observationally determined Schechter function overplot- 
ted. The simulated void galaxies well reproduce observations, and we underpredict the function in 
the "wall" sample as we our simulation volume excludes higher density regions. 




Fig. 4. — : Color magnitude diagram. On the left, the LVS (blue circles) is in reasonable agree- 
ment with the observed VGS sample (red squares). On the right, the NV (black circles) roughly 
traces the SDSS distribution for galaxies at redshifts 0.1 < z < 0.3. We see no bimodality 
i n the colors of either simulated gala xy sample such as is observed in void and field galaxies 



von Benda-Beckmann fc Muilerll2008l ). 
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Fig. 5. — : Color magnitude diagram for the NV. On the right, the sample has been split into 
red (filled circ l es) a nd blue (open circles) galaxies following the u — r = 2.2 cut identified by 



Strateva et al.1 (|200ll ). In each, the mean for the NV galaxies is over plotted in black. Also shown is 



the mean for galaxies in the void outskirts (0.5 < delta < 1.0, in green) and the mean for galaxies 
in the deepest underdensities (delta < 0.5, blue). Among the dwarf galaxies, the void galaxies are 
somewhat bluer than the galaxies in average density environments. 
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Fig. 6. — : H I mass to r-band light ratio. On the left, all simulated galaxies (b lack) are roughly an 
order of m agnitude too high com pared to the obser ved relation (in red) from IVerheiien &: Sancisi 
(|200ll ) and ISwaters et all (|2002j) (open circles) and iKreckel et all (|201ll ) (filled circles). On the 
right, the mean value for the LVS (blue), NV (black) and VS (green) shows a very slight trend with 
density within the simulated dwarf galaxies. . 
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Fig. 7. — : MHi/L r vs density for all galaxies with —15 >M r > —17, with the mean overplotted, 
reveals no clear trend with density. 




Fig. 8. — : SFR for NV galaxies (in black), with the mean overplotted for the LVS (blue), NV (black) 
and VS (green). Only at the faintest luminosities is there any difference in the three samples. 
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Fig. 9. — : S-SFR (left) and SFR/H I mass (right). The LVS (blue) and VS (green) are somewhat 
higher than the NV (blac k) for the fain t est ga laxies, but consistent with each other and with the 
void galaxies observed by lKreckel et alJ (|201ll ) (red circles). The rather large discrepancy between 
observed and simulated values in the SFR per H I mass are consistent with the large discrepancy 
in the H I masses (see Figure [6]). 
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Fig. 10. — : Mass weighted metallicity as a function of stellar mass for the LVS (blue ) with the mean 
overpl otted for LVS, NV (black) and VS (green). We reproduce the same trend as 



Tremonti et al 



(J200J) (red) for increasing metallicity with increasing stellar mass , but find no trend with de nsity. 



The observed relation, normalized to Zq at 12 + log(0/H) = 8.69 (jAllende Prieto et al.ll200ll ) with 



la error bars indicated by the shaded region, is somewhat higher as the measurements were done 
in the galaxy centers were metallicity is typically enhanced. We reproduce the mass-metallicity 
relation for the Milky Way, with Z© at M* = 5 x 1O 1O M , and the LMC, with Z©/2 at M* = 
5 x 1O 9 M . 
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Fig. 11. — : Fraction of galaxies with a given luminosity for the VC (blue dashed line), LVS (blue), 
NV (black) and VS (green). The excess of low luminosity galaxies is most apparent in the VC 
sample and relatively less dominant as a function of density , suggesting that this population is a 
result of the large scale underdensity. 




Fig. 12. — : Mean stellar ages as a function of luminosity for the LVS (blue) with the mean 
overplotted for LVS (blue solid) , VC (blue dashed) , NV (black) , and VS (green) . The low luminosity 
void galaxies are somewhat younger than the non-void galaxies. 
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Fig. 13. — : Comparison of simulated NV low-luminos ity galaxy color s with observational samples 



of dwarf galaxies from the Local Vol ume (red squares , lMakarovalll999l : iMakarova et al. 



2002, 



2009) 



and from the SDSS (blue triangles. iGeha et al.l 120061 ) . Error bars reflect uncertanties in the color 
conversion from SDSS ugzri bands to Johnson-Cousins UBVRI bands. Simulated dwarf galaxies 
are distinctly redder than observed dwarf galaxies. 
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Fig. 14. — : M r (left) and ~Mdm as a function of distance from the void center. We see no trailing 
of low luminosity galaxies into the void, however the dark matter halos do seem to segregate by 
mass. 



